Spectral function of the spiral spin state in the trestle and ladder Hubbard model 
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Eder and Ohta have found a violation of the Luttinger rule in the spectral function for the t-t' -J 
model, which was interpreted as a possible breakdown of the Tomonaga-Luttinger(TL) description 
in models where electrons can pass each other. Here we have computed the spin correlation along 
with the spectral function for the one-dimensional t-t 1 Hubbard model and two-leg Hubbard ladder. 
By varying the Hubbard U we have identified that such a phenomenon is in fact a spinless-fermion- 
like behavior of holes moving in a spiral spin configuration that has a spin correlation length of the 
system size. 

PACS numbers: 71 .1 0.-W.71 .lO.Fd 
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It is widely believed that the low-energy physics of a 
wide class of one-dimensional (ID) systems-|Can be de- 
scribed as a Tomonaga-Luttinger(TL) liquidtlu, which is 
an effective theory for electrons interacting in ID. The 
ansatz has indeed been shown to be valid for exactly solv- 
able ID models such as the Hubbard model or the super- 
symmetric t-J model, and also for some othec models nu- 
merically. However, recently, Eder and OhtaB looked into 
the spectral function in ID t-J model that has t' (t-t'-J 
model) , and made a very interesting observation that the 
density of electrons n and Fermi momentum kp are re- 
lated with kp — Ttn in a certain parameter regime, which 
is incompatible with the Luttinger relation, kp = nn/2, 
expected for a TL liquid, which they suggest to be an 
indication for a breakdown of the TL description. 

In order to clarify the origin of such an curious be- 
havior, here we study the t-t' Hubbard model with finite 
values of U. The reason we have chosen the Hubbard 
model is that the magnetic phase diagram on the n — U 
plane has been obtained for ttWrjH' Hubbard model by 
Daul and Noack (inset of Fig.2)l!m so that we can identify 
the region on which we work. For the Hubbard model we 
find the same curious behavior of the spectral function. 
We further find that the state in question is a spiral spin 
state, in which the spin correlation has a wave length of 
the system size and thus has a local ferromagnetic na- 
ture. The curious behavior of the spectral function can 
be understood by a picture in which holes can hop almost 
freely in such a spin, background as originally proposed 
by Doucot and Wend as a trial state for finite systems to 
prove the instability of Nagaoka's ferromagnetisnJa. The 
ferromagnetic-like state naturally explains the doubling 
of kp into 7m as a spinless-fermion-like behavior. The 
region over which the spiral behavior appears is indeed 
consistent with the phase diagram for the t-t 1 Hubbard 



model. In order to see if the appearance of the spiral 
state extends to other quasi-lD systems, we have also 
studied the 2-leg Hubbard ladder with U = oo, and have 
found similar features as in the t-t' Hubbard model. 
The t-t' Hubbard Hamiltonian is given by 
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in standard notations. Hereafter we set t = 1. 

First we numerically calculate, with the continued frac- 
tion expansion, the single-particle spectral function given 
by 



A ± (k,oj) 



t£I*g>. (i) 



where A + (A ) denote the electron addition (removal) 
spectrum with 7+ = c\ a , % a = ( 7j +)t, |$ G ) and £ 
the ground state and its energy, respectively. In Fig. MX 
we show the results for the case of 10 electrons on a 12 
site ring, U = 40(a) and U = 20(b) for t' = 0.2. We can 
see that for U = 20, kp = nn/2 (n = 10/12) is satisfied, 
while the equality is violated in favor of kp = 7m for 
U = 40, which would be expected for spinless fermions. 
Such a behavior persists for larger values of U. Thus 
the situation is indeed similar to tlie result for the t-t'-J 
model obtained by Eder and OhtaEl. 

At this stage, let us recall the phase diagram-pf the t-t' 
Hubbard model obtained by Daul and Noaclaffl with the 
density matrix renormalization group (DMRG) method 
in open boundary conditions. For t' = 0.2, n ~ 1 and 
large U, they found a large ferromagnetic region in the 
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phase diagram. If the periodic boundary condition(PBC) 
is assumed, on the other hand, the ground state is S = 
for any value of U at least for system sizes up to 12 sites 
as also pointed out by Daul and NoackO. We find here 
that the region in the n — U plane over which the spinlcss- 
fermion-like behavior of A(k,ui) is observed roughly co- 
incides with the ferromagnetic_region in the Daul-Noack 
phase diagram (inset of Fig. 2)3. 

Such an observation has motivated us to look into the 
spin-spin correlation function, ($G\S£Sj\<&G), to iden- 
tify the nature of the S — ground state in PBC. We 
show in Fig. || the result for 10 electrons in a 12-site ring 
with U = 20,40. We can see that for U = 40, which 
belongs to the spinless-fermion-like region, the spin-spin 
correlation indeed indicates the spiral spin state with the 
spin correlation being as large as the system size. Such 
a spiral spin-spin correlation has also been encountered 
rather ubiquitously in ID for Tasaki's mod&p, double 
exchange modelEl, the Kpndo lattice modelllil, and the 
2-band Hubbard modelEHl and akp in the 2D Hubbard 
model with U = oo for two holesEJ. 

Nature of such a spiral state has been discussed by 
Doucot and Wen(DW)d, who studied the infinite-?/ Hub- 
bard model with two holes, and found a trial state which 
gives an energy lower by inverse system size than the Na- 
gaoka's ferromagnetic state. The key idea of DW is based 
on the following intuition. Holes behave as free fermions 
for on-site interactions when the background spin state 
is ferromagnetic (or nearly so). Then we have only to 
worry about the fermion statistics, i.e., the antisymme- 
try (node) in the wave function against the exchange of 
two holes. If we impose the node upon the spin part, the 
kinetic energy of holes can be lowered. This is accom- 
plished by twisting the spin alignment in a full turn, but 
very slowly over the entire sample dimension to minimize 
the cost in energy and to maintain the ferromagnetic na- 
ture. Another way of saying is that the spiral spin tex- 
ture generates a fictitious gauge flux which absorbs the 
frustration induced by the fermion statistics of the holes. 
Indeed the flux mimics an Aharonov-Bohm of half flux 
quantum to shift the k-points by half the k-point spacing 
to let an even number of spinless fermions take a closed 
shell configuration. 

While these are conceived for the ordinary Hubbard 
model, it is intriguing to study whether the spinless- 
fermion-like behavior of A{k, u>) found here for the quasi- 
1D system can be understood quantitatively in terms of 
the DW state. Although in the original paper DW con- 
sidered the states in which the holes are dressed by spin 
waves as well, let us consider the state where holes hop 
in a rigid spin background for simplicity. The energy of 
the hole in the t-t' Hubbard model is then given as 
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where N is an integer with < N < L and ± corresponds 
to the sign of the fictitious flux. In Fig. ||a), we plot 
the energy dispersion defined by the equation (||). We 
can see that the single-particle spectrum is reproduced 
remarkably well. 

In the result for the addition spectrum an almost dis- 
persionless band of low intensity peaks is seen (at around 
E = 2 in Fig. 0(a)). If the ground state were a fully 
spin-polarized ferromagnetic state, then such a band is 
trivially expected, since we can add an opposite spin at 
any Ai-point. Remnant of such a band again suggests the 
ferromagnetic nature of the present ground state. 

A next question is whether the spiral state persists for 
more than two holes. In Fig. ||, we show the spin-spin 
correlation for 6 electrons on a 12-site lattice (n = 0.5; 
quarter filled) in PBC for U = 4, 6. We can see that 
for U = 6, the spin-spin correlation is again spiral. To 
investigate whether Doucot and Wen's picture is valid 
for such an intermediate doping, we have calculated the 
single-particle spectrum for U = 6 in Fig. [|(a), U = 4 in 
Fig. [|(b). We can see that the nature of the spectrum 
also changes between U = 4 and U — 6. The spectrum 
for the spiral state is fitted well by the energy dispersion 
defined by the eq. (||) (dashed curves in Fig. |](a)) even 
though the hole concentration is as large as rih = 0.5. 
This is surprising, since the assumption that holes are 
nearly free would be valid only for infinitesimal doping. 
For detailed energetics the sample size dependence need 
be studied. 

We next question whether the antiferromagnetic state 
is connected adiabatically to the spiral state as we in- 
crease the Hubbard U. We show the ground-state energy 
as a function of U for 6 electrons on a 12-site ring in Fig. 
§. We can clearly identify a level crossing appearing as a 
cusp around U ~ 5, which indicates that a transition oc- 
curs within the S = space from the antiferromagnetic 
phase (\ow-U side) to the spiral phase (high-?/ side). 

Finally, we move on to the 2-leg Hubbard ladder. Fer- 
romagnetic Jaehavior is also found in the ladders where U 
is infinitctirEj, so it is intriguing to study whether a spiral 
state appears in this system as well, and if so, whether 
A(k,oj) exhibits a DW-like behavior. The Hamiltonian 
is given by 
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where a(— 1, 2) labels the two legs of the ladder. We set 
t = 1, U = oo and we consider only the case of two holes. 
We have calculated the spin-spin correlation {Sf a S* a ) 
(not shown), and found that the ground state is indeed 
spiral. We have also calculated the single-particle excita- 
tion, eq.(0), where we now take 7^ = cj^ + c\ 2a (bond- 
ing states) or 7^ = c|. 1(T — c\ 2(7 (antibonding states). In 
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Fig. f|, we show the result for two holes (14 electrons 
on a 8 x 2 ladder). We fit the spectrum by a two-band 
extension of eq.(2), 

/7T\ /2Nn ?r\ 
e = -2icos(^-J cos f — — ± - J ±t. (3) 

Fig. |^ shows that the picture of two holes hopping in a 
twisted spin background is surprisingly accurate. 

To summarize, we have pointed out that the behavior 
of the spectral function found by Eder and Ohta for the t- 
t'-J model can be understood as a realization of the spiral 
spin state. As mentioned above, the fully polarized ferro- 
magnetic state is th&.ffiround state for the open-boundary 
t-t' Hubbard modeErtl This is. also the case with the 
open-boundary Hubbard ladderalil'ta. Then, in the ther- 
modynamic limit, the spiral ground state (in PBC) can 
either give way to the ferromagnetic state, remain the 
ground state, or become degenerate with the ferromag- 
netic one. This depends only on whether the cost to twist 
the spin configuration is finite or tends to zero with the 
sample size, since the energy gained by accommodating 
the electrons in a closed-shell configuration vanishes in 
the thermodynamic limit. Hence an important future 
problem is to clarify the relation between the spiral state 
and the ferromagnetic state in the thermodynamic limit 
and their spin stiffness. 

We thank R. Eder and Y. Ohta for a critical read- 
ing of the manuscript. Numerical calculation were done 
on SR2201 at the Computer Center of the University of 
Tokyo and FACOM VPP 500/40 at the Supercomputer 
center, Institute for Solid State Physics, University of 
Tokyo. 



nonmagnetic boundary is underestimated for n — > 0. For 
intermediate values of n, we can show that the ferromag- 
netic region for the anti-periodic boundary roughly coin- 
cides with that for the open boundary as will be published 
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FIG. 1. Full single particle spectral function for 10 elec- 
trons on a 12 site t-t' lattice, U=40 (a) and U=20 (b). Solid 
line denotes the electron removal spectrum, while dotted line 
the electron addition spectrum. Dashed curves in (a) indicate 
the energy dispersion defined by eq.(2). 

FIG. 2. The spin-spin correlation function for 10 elec- 
trons on a 12 site lattice, U — 20 (solid line) and 40 (dashed 
line) with the periodic boundary condition. The inset shows 
the phase boundary of the ferromagnetic phase (solid curve) 
due to Daul and Noack 2 ' 3 on which the present identifica- 
tion of the spiral phase (solid straight lines) is indicated for 
n = 10/12 (Figs. 1,2) and n = 6/12 (Figs. 3,4). 



FIG. 3. A plot similar to Fig. 2 for 6 electrons on a 12 
site t-t' lattice, U = 4 (solid line) and 6 (dashed line). 
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FIG. 4. A similar plot to Fig. 1 for 6 electrons on a 12 
site t-t' lattice, U = 6(a) and U = 4(b) 



FIG. 5. The ground state energy for 6 electrons on a 12 
site t-t' lattice as a function of 1/(7. 



FIG. 6. A similar plot to Fig. 1 for 14 electrons on a 
8x2 ladder. Dashed curves indicate the dispersion defined 
by eq.(3). 
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Fig. 1(a) 



R. Arita et al, Phys. Rev. B 



1 1 1 1 1 

(b)U=20 


i 1 






A.. 




A. 




•v 



















4.0 -2.0 0.0 2.0 4.0 

E 



Fig.1 (b) 
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Fig. 3 
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Fig.4(a) 
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Fig. 4(b) 
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Fig. 5 
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